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Abstract 

We calculate the rate of producing deep inelastic events containing an energetic isolated 
forward photon at HERA. We quantify the enhancement arising from the leading logl/x 
gluon emissions with a view to using such events to identify the underlying dynamics. 
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1. Introduction 



In the deep- inelastic small (Bjorken) x regime the gluon is the dominant partonic constituent 
of the proton, and the behaviour of the observables can be predicted in terms of the gluon 
distribution. At sufficiently small x it is necessary to resum the \ogl/x contributions which, to 
leading-order accuracy, is accomplished by the BFKL equation |lj for the gluon distribution, 
f(x,kj<), unintegrated over the gluon transverse momentum kx- As x decreases the solution 
f(x,k^) is characterized by a singular x~ x growth accompanied by a diffusion in \nk T . To 
leading logl/x accuracy, the "BFKL intercept" is A = (3as/7r) 4 In 2 for fixed as (or A ~ 0.5 
if a reasonable prescription is used to allow as to run in the BFKL equation [fj]). Sub-leading 
effects are expected to reduce the exponent. 

The singular x~ x behaviour feeds through, via the fc^-factorization theorem, into the proton 
structure function F 2 (x, Q 2 ). Thus, in principle, measurements of F 2 at HERA should be able 
to determine A. However, in practice it is difficult to distinguish the BFKL behaviour from 
that predicted by conventional DGLAP evolution, which resums as log Q 2 contributions. A 
major part of the problem is the necessity to provide non-perturbative input. Over the HERA 
regime, DGLAP evolution can mimic the steep x~ x growth by either taking a valence-like gluon 
xg(x, Qq) — ► as x — > at some sufficiently low input scale Ql [f|, or by assuming an input 
form xg(x,Ql) ~ x~ x as x — > 0, ||. Here g(x,Q 2 ) is the conventional gluon distribution 
(integrated over hr), namely 

9(x,Q 2 ) = f ^ff(x,k 2 T ). (1) 

On the other hand, in the BFKL description of F 2 , there is a diffusion in lnk T accompanying 
the translation of the input gluon f(x,k T ) to smaller x values. This leads to a significant 
contribution from the infrared region which is beyond the scope of perturbative QCD and 
which has to be included phenomenologically. 

Mueller || proposed an elegant way to overcome these ambiguities, and so sharpen up 
the theoretical description. He suggested that the study of deep-inelastic (x, Q 2 ) events which 
contain an identified energetic forward jet (xj, k 2 T ) would, in principle, allow the identification 
of the x~ x behaviour, or to be precise the (x/xj)~ x behaviour as x/xj — > 0. The process is 
illustrated in Fig. l(a)Q. Ideally we would select deep-inelastic events containing a jet whose 
transverse momentum kj? is sufficiently large to suppress the diffusion into the infrared region, 
and with longitudinal momentum fraction Xj as large as experimentally feasible so that x/xj 

2 Other properties of the final state have been proposed as possible ways to expose the mutual relationship 
between the steep power-like increase of the cross section and the lack of ordering in the transverse momenta of 
the gluon emissions along the chain. These include studies of (i) the transverse energy flow J?J , (ii) the azimuthal 
decorrelation of two jets separated by a large rapidity gap ||, |jj , and (iii) the azimuthal decorrelation of electro- 
or photo-produced dijets jl0|, O. The first tends to be masked by uncertainties due to hadronization and the 
other two by decorrelations which arise from other sources than log(l/x) resummation. 



1 



is as small as possible, typically Xj > 0.05. According to BFKL dynamics the differential 
structure function for this process has leading small x/xj behaviour 



dF 2 



<9(log l/xj)dk 2 jT 



C as{k 2 jT ) Xj 



4 / \ 

9 + g (? + <i) 




-A 



(2) 



where the parton distributions (g, q and q) are to be evaluated at (xj,k 2 T ) — where they are 
well known. The normalisation coefficient C is given in refs. fT2| , |13|, [L4], |15| . Basically the idea 



is to try to identify the BFKL behaviour from deep-inelastic scattering off a known parton at 
large k 2 T , as opposed to scattering off a proton for which we need to specify its non-perturbative 
structure at small k\. In practice, however, the clean identification and kinematic measurement 
of a forward jet, closely adjacent to the proton remnants, poses a severe experimental challenge, 
particularly for jets at the smaller values of k 2 T where the DIS + jet events are more numerous. 



Nevertheless experimental studies |16|] have been attempted with encouraging results, albeit for 
rather low values of k 2 - T . 

Here we study a related, but alternative process — namely deep-inelastic events containing 
an identified forward photon, see Fig. 1(b). As a probe of small x dynamics, the process DIS 
+ 7 has both advantages and disadvantages as compared to DIS + jet. A major advantage is 
that the measurement of a photon should be cleaner than that of a jet, and less ambiguous, 
particularly at the lower values of kr- Moreover since the q or q jet (denoted x' q in Fig. 1(b)) 
is not identified we can enlarge the data sample by including events in which its constituents 
mingle with the proton remnants. On the other hand we expect the DIS + 7 rate to be 
suppressed by a factor of order a/2ir relative to DIS + jet, though we should gain a little back 
since the photon is easier to measure than the jet and the recoiling quark jet is not identified. 
A second disadvantage is that we require isolation of the photon so as to avoid events where it 
arises from the decay of a parent ir° etc. 

We first present (in Section 2) the QCD formalism necessary to calculate the cross section 
for the DIS + forward photon process at small x. Then in Section 3 we discuss the acceptance 
criteria to be applied to the outgoing photon. In Section 4 we quantify the DIS + 7 event rate 
expected in the experiments at the HERA electron-proton collider. Finally Section 5 contains 
a brief discussion. 



2. Basic QCD formula for the DIS + 7 process 

Here we study prompt photon production in deep-inelastic scattering at small x. That is 
the process 

"7" + V — 7(z 7 > M + X > (3) 

sketched in Fig. 1(b), in which the photon is identified in the final state. Deep inelastic events 
with small x and large x 1 offer an opportunity to identify the effects of the BFKL resummation 
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of the asln(x 7 /x) contributions, which arise from the sum over the real and virtual gluon 
emissions, such as the one depicted in Fig. 1(b). In analogy to the DIS + jet process, the 
advantage of process (|3|) is that the outgoing photon acts as a trigger to select events in which 
the deep-inelastic scattering occurs off a quark in a kinematic region where its distribution, 
q(x q , k 2 T ), is known. 

The differential structure functions for this process may be written in the form 

0X^0 K^^p J ^^gT \*^Q / Q Z<yZq 

for i — T,L. The variables are indicated on Fig. 1(b) and are defined below. The region of 
integration in ([|) will be specified later. The subprocess q(k q ) +g(k g ) — > q(k' ) + ^{k 1 ) (and also 
qg — > qj) is described by the two Feynman diagrams shown in Fig. 2. It has matrix element 
squared 



\M\ 2 = 2C 2 (F) ^Sll ^ + (5) 



where C^-F) =4/3 and the invariants s and u are 

s = (ky + k' q ) 2 , u = (k q - k^f. (6) 

The fractional momenta Xi in @) are defined by the Sudakov decomposition of the particle 
4-momenta 

h = Xip' + prf + k iT (7) 

where at high energies p' and q' are to a good approximation the light-like 4-momenta associated 
with the incoming proton and virtual photon respectively 

p = p - q, q = q + xp, (8) 

2p.q 

where M is the mass of the proton. The variables Zj in (|J) and (||) are given by 

Z{ X%/Xq. (9) 

The 4-momenta of the outgoing photon and quark jet satisfy the on-mass-shell condition k 2 = 
which gives 

A = f § (10) 

Xi iq/ 

for i = 7 or the outgoing quark q' . 

The functions <&i(z, k 2 , Q 2 ) in (|J) describe the gluon chain which couples to the incoming 



virtual photon. The function <Ev corresponds to the function F in ref. [Q. To be precise, the 
factors §i/k 2 can be identified with the virtual gluon structure functions integrated over the 
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longitudinal momentum of the gluon. They can be obtained from the BFKL equation with the 
boundary condition given by the quark box (and "crossed" box) contributions. That is 



®f\z,k 2 ,Q 2 ) + 



+ ^k 2 

7T 



1 <!:' / - ilk' 2 \<5>i(z f ,k f2 ,Q 2 ) - <5>i{z f ,k 2 ,Q 2 ) <S>i(z',k 2 ,Q 2 



z Z' JO 



k' 2 



\k' 2 - k 2 \ 



\/Ak' A + k A 

(11) 



where the inhomogeneous or driving terms $^ correspond to the sum of the quark box and 
crossed-box contributions. At small z we have 



<S>T\z,k 2 ,Q 2 ) « ^\z = 0,k 2 ,Q 2 ) = ^\k 2 ,Q 2 

Q 2 



■Mi 



where 



^\k 2 ,Q 2 ) 



(12) 



2E 



1 4vr 2 



Q 4 



a s dp d 2 K (3 2 + (1 - p)' 



K 



PI 



k.(k — k) 
D\D 2 



2 E e ^«5 / dp I d 2 np 2 (i-py 

n 7T J 



l 



(13) 



with the denominators Di given by 

D 1 = k 2 + p{l - p) Q 2 



(14) 



[ K -kf + P{l-P)Q 2 . 



The BFKL equation, (|TT|), has been written for fixed QCD coupling as- In this case the 
leading small x behaviour of the solution can be given analytically. We have 



*t(z, k 2 ,Q 2 " 



9tt 2 2£e 2 a| 



(k 2 Q 2 



®l(z, k 2 ,Q 2 " 



512 ^/21C(3)/2 
2 



i ^ 



— Qp + l 



ln(l/z) 



1 + O 



\n(l/z) 



(15) 



<&t(^ 2 ,Q 2> 
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where we have omitted the Gaussian diffusion factor in \n(k 2 /Q 2 ). The Riemann zeta function 
C(3) = 1.202, and ap, the celebrated BFKL intercept, is given by 



ap 



12a, 



71 



In 2. 



(16) 



In practice we allow the coupling as to run, following the prescription of ref. ||14|| . That is 
we solve the BFKL equations for 

3a s (k 2 ^ 



H t (z,k 2 ,Q 2 " 



$i(z, k 2 , Q 2 ^ 



(17) 



Actually we solve the differential form of the equations 



dHi(z,k 2 ,Q 2 ) 3a s (k 2 ) 2 f°° dk' 2 

" rv 



HfakP^-Hfak*^) H t (z,k 2 ,Q 2 ) 



\k' 2 -k 2 \ V4fc' 4 + k A 



(18) 



<91og(l/z) 7T Jkl k' 2 

subject to the boundary conditions 

H l (z = z ,k 2 ,Q 2 ) = H?\k 2 ,Q 2 ) (19) 



where we choose Zq = 0.1 and the cut-off k 2 = 1 GeV 2 . For any small z the solutions 
Hi(z, k 2 ,Q 2 ) therefore only depend on the behaviour of Hi in the interval (z,z ). We im- 
pose the same cut, x/x q < 0.1, in the x q integration in Of course the integration is also 
subject to the kinematic constraint x q > x T We find that the slope A = ap — 1 and the overall 
normalization of the solution of the BFKL equation with running coupling as are in general 
smaller than those obtained with fixed coupling as(Q 2 )- 

We are now in a position to calculate the differential structure function dF 2 /dx 1 dk 2 T from 
dU). We compute it first using $2 obtained from solving the BFKL equation (|l8|) and then 
repeat the calculation using for $2 just the driving (quark box) term $2°^- The difference 
between the two calculations is a measure of the effect of the BFKL resummation of soft gluon 
emissions. The result for the F2 differential structure function is shown in Fig. 3 as a function 
of x 7 for three different values^ of x. We see that the difference between the BFKL and quark 
box results is dramatic once x is sufficiently small (x < 10~ 3 ) and provided x 7 ~ 0.1. Indeed 
we require x/x 7 to be small for the formalism to be valid. As explained in ref. [1^| it is the 



shape rather than the magnitude of dFi/ dx^dk 2 T which is the more reliable discriminator of 
the underlying dynamics at small x, since the normalization of the QCD predictions is subject 
to uncertainties arising mainly from the regions of low transverse momenta. 

Fig. 3 for the differential structure function for DIS + 7 events should be compared with a 



similar figure for DIS + jet events shown in Fig. 7 of ref. [^4[ . We see that the cross section is 
about a factor of 1000 lower for the photon process, as could be anticipated from the presence 
of the extra a/2iT coupling of the photon. In the next section we will impose reasonable exper- 
imental cuts on the outgoing photon, and then we will quantify the event rate which may be 
observed at HERA. 



3. Cuts to select the DIS + photon events 

We are interested in semi-inclusive deep inelastic events in which only the photon is measured 
in the final state, besides the scattered electron. On the other hand the integrals defining the 

3 The result using $2 = ^2 is independent of x (until we reach the kinematic limit /3 7 = 1), and its shape in 
x 1 reflects the quark and antiquark distributions in the proton after integration over x q . In contrast for DIS + 
jet the Xj shape of the box contribution is dominated by the gluon distribution Xjg(xj, fc| T ), that is the gluon 
is sampled "locally" . 
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differential structure functions, (f|), contain a potential singularity at s = which has to be 
regularized by a suitably chosen cut-off. A straightforward limitation of the region of integration 
to s > s , where s$ is some arbitrarily chosen value, is therefore not possible for making a direct 
comparison of (^) with experiment, since s is only determined if the recoil quark jet q' were to 
be detected as well as the photon. However, it is necessary to impose an isolation cut on the 
outgoing photon to distinguish it from background photons which arise from the decay of 7r°'s 
produced either from the proton remnants or the outgoing quark jet, q' . Suppose we impose 
an isolation angle 80 (where 80 is chosen to be around, say, 3°-10°) defining a cone around the 
photon, then at the parton level we have 8 iq > > 80 and hence a corresponding lower limit on s, 
since 

s > 2E 1 E q , (1 - cosd iq ,). (20) 

A second requirement is that our photon should be emitted in the proton hemisphere in 
the virtual photon-proton cm. frame to avoid contamination from photons radiated from the 
quark- ant iquark pair which form the quark box. Thus we require 

x 7 > /3 7 (21) 

which on using ([LC]), gives 

x 7 > ^xk 2 T /Q 2 . (22) 

The hemisphere cut, when combined with the kinematic constraint x q > x 7 , imposes an implicit 
lower limit of the x q integration in (B) 

x q > x 7 > ^xk 2 T /Q 2 . (23) 

We see that this limit is generally stronger than our BFKL requirement that x q > lOx. 

Finally there is the practical limitation that photons can only be measured if they are 
emitted at sufficiently large angles to the proton beam direction in the HERA laboratory 
frame, say 

8 1P > O . (24) 

The allowed regions of the kinematic variables (x 7 ,fc 7T ) describing the photon are shown in 
Fig. 3 for various choices of 6q. We see that large x 7 photons are only emitted at small 8 1P ; for 
a given 8 1P we can reach larger x 1 by observing photons with large k 2 T) but then the event rate 
is depleted. We also show on Fig. 3 the "forward hemisphere" boundary, (p^), for x = 6 x 10~ 4 
and Q 2 = 20 GeV 2 . 
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4. Predictions for DIS + photon production 

Here we present the numerical predictions for the cross section for the production of an 
energetic photon in the deep inelastic ep scattering taking into account various cuts discussed 
in the previous section. The relevant cross section is given by the following formula: 



da 4ira 2 



dx^dk^rpdxdC^' C^^x 



dF 2 y 2 dFr 



dX sy d k^rp 2 dX y d ^syj^ 



(25) 



We show in Fig. 5 the dependence of the integrated cross section on the variation of (a) the 
angle 9 defining the photon isolation cone and (b) the k 2 T threshold for photon detection. We 
see that the dependence on the isolation cone angle is relatively weak. 

In Fig. 6 we show the integrated cross section for prompt forward photon production as a 
function of x for three different Q 2 bins: 20-30, 30-40 and 40-50 GeV 2 respectively. We compare 
the predictions for the case where the BFKL small x resummation is incorporated with those 
where the gluon radiation is neglected. The x dependence of the cross section is driven by the 
small z behaviour of the The results show a strong enhanced increase with decreasing x 
which is characteristic of the effect of soft gluon resummation. At x ~ 10~ 4 the cross section is 
about a factor of 3.5 larger than that in which the BFKL effects are neglected. Since the impact 
factors, , are independent of x, the weak x dependence in the latter case arises mainly from 
acceptance. Finally in Fig. 7 we show the cross section in the various deep inelastic (x, Q 2 ) bins. 
The values are shown in fb and so would correspond to the number of events for an eventual 
integrated luminosity of 1 fb -1 . 

The DIS + 7 cross sections presented in Figs. 6 and 7 may be compared directly with the 
values for the DIS + jet process shown in Fig. 4 of [TS[] and Fig. 8 of JET) respectively. We see 



that there is a suppression of about a factor of 400 in going from the DIS + jet to the DIS + 
7 process. Of course the ratio depends on the precise cuts that are imposed in each case. In 
table 1 we show the effect on the DIS + 7 cross section in one of the Ax, AQ 2 bins of varying 
#0 (the minimum angle to the proton beam) and 8q (the half-angle of the isolation cone). We 
see that the event rate is less sensitive to the cone angle #0 than to the minimum angle 9q to 
the beam. 



Table 1: The DIS + 7 cross section in the bin 6 x 10~ 4 < x < 8 x 10" 4 , 20 < Q 2 < 30 GeV 2 
as calculated in Fig. 7, but for different choices of 8q and 9q. 






0o 


Aa (fb) 


5° 


3° 


46.2 


7° 


3° 


29.1 


7° 


5° 


26.2 
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5. Discussion 

The DIS + jet process is, in principle, an ideal way to probe small x dynamics, provided 
sufficiently forward jets can be measured. In practice to separate cleanly such forward jets from 
the proton remnants is a formidable challenge. Here we have studied the analogous DIS + 7 
process which has the advantage that forward photons can be more reliably measured than 
forward jets, but for which the event rate is considerably suppressed. In section 4 we quantified 
the suppression. There we predicted the DIS + 7 cross section using BFKL dynamics and 
found a characteristic rise with decreasing x, which becomes steeper the more forward are the 
detected photons. 

Our study should be regarded as an exploratory investigation of the potential usefulness 
of the process. The DIS + 7 rates which we present in Figs. 6 and 7 correspond to photons 
produced (i) at more than 5° to the proton beam direction (9 1P > 60 = 5°) in the HERA frame, 

(ii) in the centre of an isolation cone of half-angle 3° {Q iq > > 9 = 3°) in the HERA frame, 

(iii) in the proton hemisphere in the 7*p cm. frame (to avoid contamination with photons 
radiated from the quark box), and (iv) with k^ T > 5 GeV 2 . Our work is at the parton level 
and so in practice the isolation criteria will need further study. In particular simulations of 
the fragmentation of the outgoing (q') jet should be performed so as to be able to choose 
the optimum isolation criteria for the photon. Recall that isolation is required to suppress 
(background) photons from 7r° decays. 

There is a second reason why it is important to know the rate of DIS + forward photon 
events. We have stressed the experimental difficulty of obtaining a clean sample of DIS events 
containing a measured forward jet close to the proton remnants. An alternative possibility to 
expose \og(l/x) resummation effects is to use the improved knowledge of the fragmentation 
functions to identify the forward jet through the measurement of a single energetic decay prod- 
uct. As it happens the n° is the hadron which can be identified in the most forward direction 
in the detectors at HERA. For this alternative signal of small x dynamics, our process is then 
the background. 
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Figure Captions 



Fig. 1 Diagrammatic representation of (a) a deep-inelastic + forward jet event, and (b) a deep- 
inelastic (x, Q 2 ) + forward identified photon (x 7 , k^r) event. 

Fig. 2 The Feynman diagrams describing the gq — > 79 subprocess embodied in the DIS + 7 
diagram shown in Fig. 1(b). 

Fig. 3 The F 2 differential structure function, @, for deep-inelastic (x,Q 2 ) events accompanied 
by a measured forward photon (x 7 , k 2 T ) as a function of x 1 for different values of x, x — 
10 -4 , 10~ 3 and 10~ 2 , and for Q 2 = 5 GeV 2 . The continuous curves correspond to inputting 
into (|4]) the solution $ 2 of the BFKL equation ([18]), whereas the dashed curve is calculated 
using for $2 simply the driving term $2 °f (P^D- Here \A4 \ 2 is regulated by requiring s > 1 
GeV 2 . Plots (a) and (b) correspond to photons with k 2 T = 5 and 10 GeV 2 respectively. 

Fig. 4 The curves give the upper boundary of the allowed regions of the photon kinematic vari- 
ables (xj, k 2 T ) for deep-inelastic + photon events with x = 6 x 10~ 4 and Q 2 = 20 GeV 2 
for various choices of the acceptance angle 9 of (^4j). The photon angle 9 to the proton 
direction in the HERA (30 x 820 GeV) laboratory frame is not uniquely specified by 
(x, Q 2 ; Xy, k 2 T ). Varying the remaining azimuthal angle transforms the line of constant 
8 1P into narrow bands in the x 7 , k 2 T plane. The continuous lines that are shown are 
obtained by averaging over the azimuthal degree of freedom. The lines are insensitive to 
reasonable variations of x and Q 2 . Also shown (by a dashed line) is the lower boundary 
given by the hemisphere cut ( [2"2"[ ) for x = 6 x 10~ 4 and Q 2 = 20 GeV 2 . 

Fig. 5 The dependence of the DIS + 7 cross section integrated over the photon variables to 
variation of (a) the angle 8 defining the isolation cone around the photon {9 1 j > 8 ) 
and (b) the threshold for photon detection (k 2 T > k 2 h ). The results are for the x, Q 2 bin 
defined by 6 x 10~ 4 < x < 8 x 10~ 4 and 20 < Q 2 < 30 GeV 2 . We impose the hemisphere 
cut (p3|) and take 9 = 5° in (fj4|). (The lack of smoothness of the curves simply reflects 
the errors on the six- fold numerical integration). 

Fig. 6 The cross section, (a) in pb, for deep inelastic + photon events integrated over Ax = 
2 x 10~ 4 , AQ 2 = 10 GeV 2 bins which are accessible at HERA, and integrated over the 
region 1P > 5°, k 2 T > 5 GeV 2 , but subject to ( [23] ) and an isolation cut of 6q = 3°. The 
x dependence is shown for three different AQ 2 bins, namely (20,30), (30,40) and (40,50) 
GeV 2 . The (a) values are plotted at the central x value in each Ax bin and joined by 
straight lines. The continuous curves show (a) calculated with determined from the 
BFKL equation, whereas the dashed curves are obtained just from the driving terms 
i.e. from the quark box. For clarity a vertical line links the pair of curves belonging to 
the same AQ 2 bin. 
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7 The cross section (a) in fb for deep inelastic + photon events in various (Ax, AQ 2 ) bins 
which are accessible at HERA, and integrated over the region Q yp > 5°, k^ T > 5 GeV 2 , but 
subject to ( [25D and an isolation cut of 6> = 3°. The number in brackets is the cross section 
calculated with just the quark box approximation (<E>j = The difference between 

the two numbers is therefore the enhancement due to the BFKL soft gluon resummation. 
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(a) DIS + jet (b) DIS + y 
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